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ABSTRACT 

The stars in a Globular cluster (GC) have always been considered coeval and of the same metallicity Recently, this assumption has 
been challenged on the basis of spectroscopic and photometric observations, which show the existence of various generations of 
stars in GCs, differing in the abundances of products of H-buming at high temperatures. The main final product of this burning is 
He. It is then important to study the connections between stars properties and He content. We consider here the about 1400 stars on 
the Red Giant Branch (RGB) observed with FLAMES@VLT in 19 Galactic Globular Clusters (GCs) in the course of our project 
on the Na-O anticorrelation. Stars with different He are expected to have different temperatures (i.e., different colours), slightly 
different metallicities [Fe/H], and different luminosity levels of the RGB bump. All these differences are small, but our study has the 
necessary precision, good statistics, and homogeneity to detect them. Besides considering the observed colours and the temperatures 
and metallicities determined in our survey, we computed suitable sets of stellar models - fully consistent with those present in the 
BaSTl archive - for various assumptions about the initial helium content. We find that differences in observable quantities that can 
be attributed to variations in He content are generally detectable between stars of the Primordial (P, first-generation) and Extreme (E, 
second-generation) populations, but not between the Primordial and Intermediate ones (1). The only exception, where differences are 
significant also between P and 1 populations, is the cluster NGC 2808, where three populations are clearly separated also on the Main 
Sequence and possibly on the Horizontal Branch. The average enhancement in the He mass fraction Y between P and E stars is about 
0.05-0.11, depending on the assumptions. The differences in Y, for NGC 2808 alone, are of about 0.11-0.14 between P and 1 stars, 
and about 0.15-0.19 between P and E stars, again depending on the assumptions. When we consider the RGB bump luminosity of first 
and second-generation stars we find different levels; the implied Y difference is more difficult to quantify, but is in agreement with the 
other determinations. 

Key words. Stars: abundances - Stars: atmospheres - Stars: Population II - Galaxy: globular clusters - Galaxy: globular clusters: 
individual: NGC 104 (47 Tuc), NGC 288, NGC 1904 (M 79), NGC 2808, NGC 3201, NGC 4590 (M 68), NGC 5904 (M 5), NGC 6121 
(M 4), NGC 6171 (M 107), NGC 6218 (M 12), NGC 6254 (M 10), NGC 6388, NGC 6397, NGC 6441, NGC 6752, NGC 6809 (M 55), 
NGC 6838 (M 71), NGC 7078 (M 15), NGC 7099 (M 30) 



1. Introduction 

. 5__, ' Recent progresses have indicated that large star-to-star variations 

S^ exist in the He content among stars in globular clusters (GCs), 

;_( in spite of difficul ties in d irect deriva tion of He abund ances (see 

_C3 ^ however .Villanova et al.L 12009: .Moehler et al.L l2007h . The first 

clear evidence has been obt ained from observa tions of multiple 

main sequences in w Cen ( BedinetaU 20041) and NGC 2808 

JD^^ntona et al, 2005; Pi otto et al.L l2007h . In co Cen, stars on 

the bluest main seq uence (MS) are m ore metal-rich than those 

of the reddest one dPiotto et al.L l2005t) . a fact that can only be 

explained by an He content higher by about 0.10 to 0.15 in the 

mass fraction Y (as originally suggested by NoiTis, '2004") with 

respect to a canonical (primordial) He content of about 0.245 

commonly adopted for the bulk of the stellar population in this 

cluster The small scatter found in NGC 2808 among red giant 



branch stars also implies a larger He content for the bluer MS, by 
roughly the same amoun t. Two other clus t ers ha ve been found to 
have wide MSs: 47 Tuc dAnderson et al.L 12009) and NGC 6752 
dMilone et al.Ll2010l who also suggest a possible split). Q 

A star-to-star spread in the He abundance may explain many 
aspects of the horizontal br anches (HB) of GCs, as extensively 
discussed in Gratton e t aP (l2010h . but we are aware that also 
the case fo r a null s pread of He in GCs is maintained in liter- 
ature (e.g., ICatelanl l2009l) . Internal variations in He are likely 
connected with other chemical signatures found in GCs, related 



* Based on observations collected at ESO telescopes under pro- 
grammes 072.D-507 and 073.D-0211 



' T he split subgi ant branches (SGB e.g., NGC 6388: iMoretti et aP 
120091 ; NGC 1851: iMilone et al.l 1200 8^ seem more related to differ- 
ences in age or CNO content. A spread in age is probably the cause 
of the spread/split MS turn-off's found in many intermediate-age 
Magellanic Clouds GCs (Mackevetal. 2008; Milone et al. 2009 in 
LMC; Glatt et al. 2008 in SMC) even if for them -but notably not for the 
old Galactic GCs- there is an alternative physical interpretation based 
on the observational effects associated to the o ccurrence of stellar rota- 
tion as discussed bv lBastian & de MinkI i2009h . 
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to H burning occurring at hig h temperature. This has been first 
proposed in a modern way bv lD'Antona et al] (l2002 l). who no- 
ticed that the smaller masses of He-rich stars may naturally ex- 
plain the correlation between the presence of stars at the blue 
extreme of the horizontal br anch, and of stars extremely de- 
pleted in O in various clusters. TLee et alj (12005) suggested that a 
super-He rich population could explain the extreme HB stars in 
GCs, specifically presenting the cases of co Cen and NGC 2808 
(where they predicted a spread or a split of the MS, just before 
it was actually found). iKavirai et al.l (l2007 i) suggested that the 
far-UV properties of GCs in the Virgo elliptical M87 could be 
explained by a su per-He rich population similar to the one in co 
Cen. lD' Antona& Caloil (|2008) tried to explain the whole dis- 
tribution of stars along the HB of various clusters (including 
NGC 2808, M 13, and M 3) as d ue to star- to-star variations in the 
He co ntent (but see a lso Castella ni et alJ l2005). In a companion 
paper dGratton et al.L 120101) . we make a more systematic com- 
parison of the extremes of the distribution of stars along the HB 
and of the Na-O anticorrelation. We find that median HB colours 
are quite well reproduced by a combination of metallicity, age, 
and a simple mass loss proportional to metallicity. However, the 
low-mass extreme of the HB is clearly cor related with the exten- 
sion of the Na-O anticorrelation (see also lCarretta et aUl2009al) 
which, in turn, seems to be linked with the cluster absolute mag- 
nitude, a proxy for cluster mass. 

These correlations suggest that the multiple generation sce- 
nario, needed to explain the Na-O anticorrelation]3 is able to ex- 
plain many peculiarities of GCs. As we have seen, a crucial da- 
tum in this scenario is the He content of second-generation stars. 
While the most spectacular evidence of a variation of the He con- 
tent comes from the HB and the MS, also the red giant branch 
(RGB) stars are very useful because chemical composition (but 
-alas- not He content) can be determined easily enough for these 
stars using spectroscopy, unlike faint MS stars. This allows us 
to cross correlate the abundances of various elements, provid- 
ing important constraints on the properties of the He producer, 
which have not been clearly identified yet (two main possible 
mechanisms have been proposed, intermediate/massive asymp- 
totic gia nt stars -AGB- and fast rota ting massive stars -FRMS- 
see e.g.. iD'Antona & Ventural l2007l and iDecressin et al.ll2007l 
respectively). Various properties of RGB stars are affected by 
the variable He abundances: 

1. The temperature (at a given luminosity) of He-rich RGB 
stars is shghtly warme r than that of He-poor ones (see e.g., 
ID'Antona et al.L l2002h . The difference is not large (only a 
few tens of K), and very accurate relative determinations are 
required. However, both colours and line excitation can be 
used to show that this really occurs. 

2. The gravity (again at a given luminosity) is also different, due 
to both the smaller mass and higher temperature of He-rich 
stars. He-rich stars have slightly lower gravity, but this effect 
is very small, and difficult to show up, because of several 
complications (for instance, the larger molecular weight and 
reduced continuum opacity simulate a higher gravity, offset- 
ting most of the effect expected on e.g., ionisation equilib- 
rium) so we will not consider it here. 



^ It is generally assumed that first-generation stars have chemical 
composition similar to field stars of similar age and metallicity (i.e., 
high O and Mg, low Na and Al, etc.), while second-generation stars, 
formed from material polluted by matter lost from the primordial popu- 
lation, show the abundance patterns (e.g., Na-O anticorrelation) that are 
unique of GCs. 



3. The [Fe/H] value is different, because of the different H con- 
tent (broadly speaking, if Y increases, X decreases, hence 
[Fe/H] increases; things are of course more complicated, see 
Sec. 12.3b : this difference is small but n ot negligible. We al- 
ready noticed this effect in NGC 2808 JCarretta et all l2006l 
Paper I). 

4. Finally, the luminosity of the RGB bump is also different: 
the R GB bump is more luminous for He-rich stars ( Ibea 
1968; ISaIaris et al.[ 120061) . even if also in this case the situa- 
tion is more complex (see Sect. [3]). We showed that possible 
evidence for this effect can be traced in the run of Na abun- 
dances along the RGB of NGC 6218 (=M 12) andNGC 6752 
(Carrettaetal.. 2007b. Paper IV). 

In this paper, we intend to look for evidence of a variable 
He content in the sample of more than 1,400 stars along the 
RGB of 19 GCs for which we acquired FLAMES/GIRAFFE 
spe ctra in the last few y ears (presented in a series of papers: 
see ICarretta et aI.L l2009alfbL -Paper VII and VIII- and references 
therein). Several features of this data set make it very promis- 
ing for the present purpose: the sample is very large; it has been 
analysed using a very homogenous technique; star-to-star errors 
in effective temperatures, usually the main source of errors in 
this analysis, were reduced to very small values by a carefully 
tailored technique. In addition, we have determined Na and O 
abundances for almost all these stars, allowing to classify them 
into three different groups according to Na enhancement and 
O depletion (Paper VII). We called these groups P=primordial, 
that is stars with composition similar to that of field stars, and 
likely belonging to the primordial population; I=intermediate, 
and E=extreme, that are stars with different degrees of O deple- 
tion and Na (and likely also He) enhancement. We recall here 
that we conservatively applied this definition only to stars with 
both O and Na measured, thus reducing our sample to about 960 
objects, the ones that will be used in the present paper On the 
other hand, the way we analysed stars should be carefully con- 
sidered when we try to extract information from our data set. 
The most significant assumptions we made are that stars of the 
same luminosity have the same effective temperature and surface 
gravity. Whenever needed, we will consider explicitly the con- 
sequences of these assumptions, and we will correct our results 
accordingly. 

Beside observations, determination of the impact of different 
He abundances requires appropriate modellin g. To this purpose 
we used a set of BaSTI evolutionary models (Pietri nferni et al.L 
2004, 2006, 2009) purposely computed for this paper employing 
two sets of heavy elements mixtures (see Sect.|2]and[3]for details 
and relations). 

The structure of this paper is as follows: in Sec.|2]we present 
the observed differences in colour, temperature, and metallicity 
and derive the implied differences in He content; in Sec. [3] we 
compare the RGB bumps of first and second-generation stars, 
deducing a different He abundance also in this case. Finally, we 
summarise and discuss our results in Sec. |4] 



2. Difference in He from RGB stars 

We may find different He values for the three populations us- 
ing various indicators. To actually compute the implied Y (as in 
the definition Xh-Y-i-Z=1 for the chemical composition of stel- 
lar models on mass fraction), we used relations based on the 
BaSTI stellar evolutionary models. We have computed an ex- 
tended set of low-mass stellar models for various assumptions 
about the initial He content. All the stellar model predictions 
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Fig. 1. Mean differences in {V - K) colours between the I and P 
components in the 19 GCs as a function of the metallicity dif- 
ference (I-P). The two are well correlated, as indicated by the 
Pearson and Spearman coefficients. 



adopted in the present analysis have been obtained exactly in 
the sam e physical fram ework used for the BaSTJj stellar model 
library dPietrinfemi et a l.. 2004, 2006). Thus, they are fully con- 
sistent with the predictions corresponding to a 'canonical' as- 
sumption about the initial He abundance (0.245 < F < 0.27). In 
particular, we used two assumptions for the heavy element mix- 
tures (others may be valid, but these represent two extremes): 

(fl) one where the second-generation I, E stars have the same 
heavy elements distribution of the first-g eneration P stars, 
both a-enhanced (for more details see iPietrinferni et al.L 
|2006) and without peculiarities in the distribution of C, N, 
O, Na (a simplified approach); 

{b) a second one where the I, E stars have a peculiar distribu- 
tion of metals (i.e., showing a signature of the Na-0 an- 
ticorrelation) and different from the standard one of the P 
stars. In particular, we adopted the "CNONa-extreme" chem- 
ical composition accounted for by Pietrinferni et al. (2009), 
which has the C+N-i-O sum enhanced by a factor of two. 

It is worth mentioning that, for a fixed iron content [Fe/H], 
the global metallicity Z is higher for case b than for case a 
mainly because the sum of (C-(-N-i-0) is enhanced with respect 
the 'normal' a-enhanced mixture. This has some implications: 
for the same He, the temperature of a CNONa-extreme track is 
lower than the corresponding normal a-enhanced and its RGB 
bump is fainter (see Sect. 3). This also means that we cannot 
give a simple and straightforward interpretation of differences 
e.g., in colour or RGB bump brightness as differences in He be- 
cause the result also depends on the mixture of heavy elements 
we attribute to the second-generation stars. 



^ The whole set of stellar models used in present work, as well as 
additional predictions for He-enhanced models can be retrieved from 
the URL site: http : //www . oa-teramo . inaf.it/BASTI i 
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Fig. 2. Individual differences I-P of metallicity (upper panel) 
and colours (lower panel) versus [Fe/H], with Spearman and 
Pearson coefficients indicated (in the lower panel, values in 
parenthesis are obtained excluding NGC 6388 and NGC 6441). 



2.1. A y from difference in colours 

For all stars in our 19 GCs we have Johnson V magnitudes. 
The photometr y was obtained as described i n P apers I to VIII 
of our se riesgf 'Carretta et alj (l2006i l2007allblQ |2009a"b') and 
iGratton et al. ( 200 6, 2007), where the interested reader can find 
references to the original papers or description for the unpub- 
lished data. We have retrieved K magnitudes from the 2MASS|_| 
Point Source Catalog dSkrutskiee t al., 2006) for almost all stars. 
Distance moduli and colour excesses E{B - V) were taken from 
the on-line version of the GC catalog bv Harris (1996) . 

In our series of papers devoted to the study of Na and O 
abundances for RGB stars, we derived the temperatures of the 
stars using the Alonso et al. (1999) relations and the photome- 
try. In particular, to decrease the internal error, we derived the 
temperatures not directly from the {V - K) colours, but from re- 
lations of temperature with the V or A" magnitude (see Paper 
II to VIII; the exception is NGC 2808, the first cluster anal- 
ysed, for which we used the classical relation between colour 



^ In particular: in Paper I for NGC 2808, Paper II for NGC 6752, 
Papers III, V for NGC 6441, Paper IV for NGC 6218, Paper VI for 
NGC 6838, Paper VII for NGC 104, NGC 288, NGC 1904, NGC 3201, 
NGC 4590, NGC 5904, NGC 6121, NGC 6171, NGC 6254, NGC 6397, 
NGC 6809, NGC 6838, NGC 7078, and NGC 7099 

^ The Two Micron All Sky Survey is a joint project of the 
University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 
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Table 1. Difference in colour and metallicity between different populations in the 19 GCs. 
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mean 
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mean 
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mean 


rms 




NGC 104 


-0.768 


24 


0.033 


0.061 


52 


-0.010 


0.038 


-0.760 


0.041 


-0.748 


0.039 


0.012 


NGC 288 


-1.305 


20 


-0.004 


0.053 


41 


-0.013 


0.048 


-1.230 


0.046 


-1.223 


0.052 


0.007 


NGC1904 


-1.579 


17 


-0.010 


0.050 


16 


-0.008 


0.063 


-1.549 


0.042 


-1.536 


0.042 


0.013 


NGC2808 


-1.151 


20 


0.001 


0.059 


11 


-0.046 


0.027 


-1.167 


0.033 


-1.097 


0.056 


0.070 


NGC3201 


-1.512 


31 


-0.017 


0.063 


54 


0.003 


0.063 


-1.468 


0.061 


-1.494 


0.053 


-0.026 


NGC4590 


-2.265 


18 


-0.014 


0.043 


28 


0.002 


0.040 


-2.215 


0.065 


-2.236 


0.052 


-0.021 


NGC5904 


-1.340 


21 


0.002 


0.046 


53 


-0.006 


0.043 


-1.342 


0.042 


-1.343 


0.034 


-0.001 


NGC6121 


-1.168 


23 


0.179 


0.083 


53 


0.008 


0.070 


-1.310 


0.049 


-1.302 


0.045 


0.008 


NGC6171 


-1.033 


10 


-0.013 


0.058 


16 


0.010 


0.087 


-1.072 


0.053 


-1.054 


0.072 


0.018 


NGC6218 


-1.330 


18 


0.000 


0.047 


43 


-0.024 


0.053 


-1.302 


0.028 


-1.302 


0.041 


0.000 


NGC6254 


-1.575 


31 


-0.011 


0.058 


48 


-0.009 


0.064 


-1.550 


0.072 


-1.552 


0.057 


-0.002 


NGC6388 


-0.441 


12 


0.007 


0.097 


6 


-0.035 


0.062 


-0.427 


0.115 


-0.389 


0.061 


0.038 


NGC6397 


-1.988 


















-2.006 


0.067 




NGC6441 


-0.430 


3 


-0.062 


0.084 


3 


0.047 


0.043 


-0.284 




-0.311 




-0.027 


NGC6752 


-1.555 


21 


-0.037 


0.053 


57 


0.027 


0.055 


-1.549 


0.055 


-1.563 


0.053 


-0.014 


NGC6809 


-1.934 


16 


-0.035 


0.062 


55 


0.017 


0.053 


-1.946 


0.042 


-1.955 


0.057 


-0.009 


NGC6838 


-0.832 


11 


-0.001 


0.042 


19 


-0.005 


0.041 


-0.803 


0.044 


-0.808 


0.038 


-0.005 


NGC7078 


-2.320 


12 


-0.029 


0.032 


18 


0.003 


0.069 


-2.312 


0.066 


-2.320 


0.050 


-0.008 


NGC7099 


-2.344 


11 


-0.039 


0.048 


13 


0.023 


0.059 


-2.313 


0.046 


-2.340 


0.053 


-0.027 



and temperature). In the present paper we used the unreddened 
(V - K) colours. The second step provided the necessary elim- 
ination of offsets between the photometries: in each cluster we 
selected stars belonging to the P component, computed the aver- 
age offset in colour between (V - K) observed and derived from 
the colour-temperature relations and used it to define the zero 
point for the P stars (so that they have zero colour by definition). 

Only after this normalisation we computed the offsets in 
colours between the three populations. We have 319, 587, and 
39 stars in the P, I, E components, respectively. We only con- 
sidered stars fainter than Mk - -3.5 because for brighter stars 
the different sequences for different He are too close and become 
indistinguishablqj- 

Results for the 19 GCs are given in Table [l] where we indi- 
cate the cluster, the metallicity (Col. 2, taken from Paper VIII), 
the number of stars in the P component (Col. 3) and their average 
colour (y - K)p and r.m.s. (Cols. 4, 5), the number of stars in the 
I component (Col. 6) and their average offset with respect to the 
P stars, with rm.s. (Cols. 7, 8); Cols 9 to 12 show the average 
metallicities of the P and I components, with their r.m.s., and 
Col. 13 shows the difference in [Fe/H] (see next Section). We 
did not evaluate the corresponding individual cluster values for 
the E component since the number of E stars is very small and 
they would not be significant. Figure [1] shows that differences in 
colours and in metallicity are well anticorrelated as indicated by 
the Pearson and Spearman rank coefficients; this is significant at 
better than the 99.9% level. 

To see whether there is a difference in colour among the dif- 
ferent populations we summed up the P, I, and E components in 



* This has also the effect of minimising the possible contamination of 
AGB stars, since at that magnitude level the separation between RGB 
and AGB is clear. As an example, in the c ase of NGC 6752 , M^ = -3.5 
means V~12.4 fsee lCarretta et alil2007ah . From Fig. 1 of that paper, the 
separation in colour (B - V) between the two sequences is then about 
0.1 mag, i.e., much more than the photometric error. 



all clusters and obtained the following weighted meano 
(V - K)i -iV- K)p = -0.003 + 0.004 {r.m.s. = 0.060) 
(y - K)e - (y - K)p = -0.037 + 0.010 {r.m.s. = 0.056) 

We see that A(y - K) is of the same order of the eiTor, hence not 
significant, for the P and I stars, while it is about 4 times the error 
for the P and E stars. We conclude that the extreme populations 
is significantly bluer than the primordial one. This A(y - K) cor- 
responds to a AY of about 0.05 or 0.08, according to the models 
cited, for case a) and b), respectively. The differences in colour 
and in metallicity (see Sec. 2.3) are summarised in Table|2l while 
the corresponding AY values are shown in Table [3] for simplic- 
ity. 

A complication may arise from the presence of CO bands 
in the K filter They are stronger in P stars than in 0-depleted 
I and E ones, and may depress the flux in the K filter by about 
1-2%, making the P stars to appear bluer, hence decreasing the 
difference in {V—K) among populations. From simulations based 
on the CO index by ICohen et alJ (Il978h we have estimated that 
the (y-/r) of P stars should be about 0.013 mag redder than they 
appear. That means that the corrected A(y - K) between E and P 
is about -0.05; this corresponds to AY of about 0.06 or 0.10, for 
the two mixtures. 

This is a global result, averaged over the 19 GCs; it is inter- 
esting to consider NGC 2808 separately, the only cluster in our 
sample for which large differenc es in He have bee n deduced on 
the basis of the multiple MSs tPiotto et all l2007h . We find for 
NGC 2808: 

(y - K)i - (y - K)p = -0.036 + 0.015, 

(y - K)e - (y - K)p = -0.044 + 0.017 

For this cluster both the I-P difference and the E-P one are sig- 
nificant at about the 2.5cr level. These A(y - K) imply a AY 



' The values would become -0.003 ± 0.004 and -0.042 ± 0.008 had 
we eliminated NGC 3201, a cluster with di fferential redde ning (which 
we corrected for anyway, using the maps bv lvon Braun & M ateo 20010. 
Its potential influence on the global properties of our sample is larger 
than the one by e.g., NGC 6388 or NGC 6441 because we have a much 
larger number of stars. Even so, the results are indistinguishable. 
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Table 2. Summary of differences in colour and metallicity be- 
tween populations. 



Sample 


HV-K),p 


MV-Khf. 


A[Fe/H],p 


A[Fe/H]EP 


all 


-0.003 


-0.037 


0.000 


0.027 


all(+corr)' 


- 


-0.050 


- 


- 


NGC 2808 


-0.036 


-0.044 


0.046 


0.084 


NGC 2808(+corr)' 


-0.049 


-0.057 


- 


- 



Notes. With correction for the CO bands. 



of about 0.05 (case a) or 0.08 (case b), and 0.08 (0.09) be- 
tween the P and I or E populations, respectively. Taking again 
into account the effect of CO bands, these differences become 
A(y - K)i^p = -0.049 and A(V - K)e-p = -0.057 which trans- 
late into AY/_p = 0.06 (0.10) and AY e-p = 0.07 (0.11), respec- 
tivelyQ 

2.2. AY from difference in dO 

In alternative to the use of temperatures derived from colours, we 
may also use information derived from excitation to find (small) 
temperature differences between the three populationsQ This is 
possible because we did not derive T^f/s from the excitation 
equilibrium (i.e. the spectroscopic route) but from the photome- 
try. From the excitation equilibrium we derived the d6 values for 
all the programme stars; recalling that d9 > Q means that the star 
is warmer than expected from its luminosity, we found the fol- 
lowing weighted average values: d6E - d0p - +0.0092 + 0.0040 
using the 39 E stars and the 207 P stars with more than 25 lines 
measured. The result improves to d6E - dOp - +0.0094 ± 0.0037 
if we relax the request on the number of lines and use the 61 E 
stars and the 315 P stars with more than 15 lines measured. This 
difference is significant at 2.5cr and corresponds, at 4500 K, the 
mean temperature of the stars considered, to a difference in tem- 
perature of 32 ± 14 K, in the sense of E stars being warmer than P 
ones. This implies a AY of about 0.06 between E and P stars (an 
average between results for the cases a) and b) discussed above). 
Since in this case we have averaged stars of all magnitudes and 
the sequences tend to converge at brighter luminosity, the aver- 
age difference is smaller than the one implied by colours. 

Once again, the difference between the P and I populations 
is not significant; the corresponding differences are d9i - dOp - 
-0.0037 ±0.0017 and £/0/-£/0p = -0.0018 ±0.0018 considering 
stars with at least 25 or 15 lines measured, respectively. 

2.3. AY from difference in [Fe/H] 

Grossly speaking, given the definition X+Y+Z=l, for a fixed 
global metallicity Z, an increase in He (Y) has to be associated 
to a decrease of the H abundance, so that [Fe/H] increases. Of 
course, the issue is more complex, since a change in the H to 
He ratio has effects on the star structure, which are taken into 
account by the evolutionary models used i n the present paper. 

We noticed this feature in NGC 2808 (ICarretta et al.L l2006t) 
when, dividing stars in a way different from the one adopted 



^ The small difference in AY£_p between the average of 19 GCs and 
NGC 2808 alone is due to the fact that the latter has a strong contribu- 
tion to the global E population. 



The definition is: 6, 



5040/r„c. Since r„r - 0.86 x T,ff, we 



have dgjc = 5S60 /Tgff. Differentiating the equation we have the relation 



between dO and dr, 



'eff 
'<■//■ 



here, we found slightly different (and increasing) values of 
[Fe/H] for "0-normal", intermediate, and very O-poor stars (i.e., 
with canonical, intermediate, and enhanced He in the usual inter- 
pretation of Na-O anti-correlation). We explore here the whole 
sample of our programme clusters. 

In this case we have to proceed in a slightly different way. As 
mentioned in Sect. 2. 1, in our series of papers (Paper II to Paper 
VII, but notably not in Paper I on NGC 2808) we derived the 
temperatures used to obtain [Fe/H] and other abundance ratios 
from a calibration between magnitude and colour-based tem- 
perature. This reduced the effect of errors (since magnitudes 
are more reliably measured than colours) and allowed to have 
smaller internal eiTors, but also effectively collapsed all possible 
differences in metallicity due to different He, that would have 
shown as different colours, hence temperatures. Before compar- 
ing [Fe/H] values of the different populations we have then ap- 
plied a correction to transform the metallicities to the values they 
would have had, had we used directly the colour-based T^j/s. 

When we apply the same procedure adopted for colours, i.e. 
normalising to the P population of each cluster (for a total of 320 
stars), we obtain the following weighted means: 

[Fe/H]i - [Fe/H]p = 0.000 ± 0.003 (r.m.s = 0.051) 

[Fe/H]E - [Fe/H]p = 0.027 ± 0.010 (r.m.s = 0.059) 

with 587 stars in the I and 38 stars in the E component, respec- 
tively. If we eliminate NGC 3201 because of its differential red- 
dening, things do not change much: these two numbers become 
0.002 ± 0.004 (288 stars, rm.s=0.055) and 0.032 ± 0.010 (32 
stars, rm.s. =0.059), respectively. Only the difference between E 
and P populations is significant and, for the two cases, it implies 
a AY ^ 0.05 (or 0.06, without NGC 3201) if the whole differ- 
ence is due to the denominator in [Fe/H], i.e. if the distribution 
of heavy elements is the same in P and E stars (case a). The 
situation is more complicated in case b, where also Z changes; 
we have run models with different assumptions for the change in 
Z (from zero to 0.002); a A[Fe/H]=0.03 implies for an average 
difference in Z of 0.0013, a AY ^ 0.11. 

Also in this case we may separate NGC 2808 and we obtain 

[Fe/H]i - [Fe/H]p = 0.046 ± 0.020 (r.m.s = 0.066) 

[Fe/H]E - [Fe/H]p = 0.084 ± 0.018 (r.m.s = 0.052) 

Both values are significant and imply AYj^p = 0.11 and 
AY E-p - 0.15 (case a). For case b, these values become, again 
assuming the same difference in Z as above, AY/_p = 0.14 and 
AY E-p =0.19 



2.4. Differences among fhe tiiree populations 

Even if this is part of our global study of the anticorrelations 
in GCs, we have to remember that AY is not a direct synonym 
of anticorrelation and that the same extension of the anticorrela- 
tion s (as measured, e.g. from the Interquartile range IQR[0/Na], 
see ICarretta et alJ l2006l l2007dh . may be reached for different 
values of AY. We have an indication of that from Fig. |2l in the 
lower panel we plot the difference in metallicity and colour be- 
tween I and P stars as a function of cluster metallicity, while 
in the upper panel we do the same for differences in metallic- 
ity. The individual values have large eiTor (and this is why we 
considered only the averages over the whole sample); however, 
there is a hint that the effects of differences in He content (as 
inferred from the difference in metallicity or colour) are less evi- 
dent for the metal-poor clusters. Only in the metal-richer clusters 
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Table 3. Summary of differences in Y between first and second-generation stars. 



Sample 



AYpi(V-A:) AYpe{V-K) AYpi[Fe/H] AYpE[Fe/H] AYpi(V-A:) AYpe(V - K) AYn[Fe/H] AYpe[F<?///] 



Case a 



Case b 



all 




0.049 


all(+corr.)' 




0.064 


NGC 2808 


0.047 


0.057 


NGC 2808(+corr.)' 


0.063 


0.072 



0.109 



0.046 
0.146 



0.075 
0.098 



0.077 
0.100 
0.089 
0.113 



0.14 



inr 

0.19 



Notes. With correction for the CO bands. 
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Fig. 3. Cumulative distributions for difference in colour (upper 
panel) and metallicity (lower panel) for individual stars in the 
three populations. 



the I stars are bluer and metal-richer than the P ones. This could 
suggest that, for the same production of proton-capture elements 
(e.g., destruction of O and production of Na), more He is pro- 
duced in metal-richer clusters. In other words, in the metal-poor 
GCs a smaller AY is required to produce an effect on the shape 



of the Na-0 anticorrelation, i.e., on the classification of a star in 
the I or even E population. 

While the effect discussed above needs corroboration, a solid 
conclusion is that, with the exception of NGC 2808, the differ- 
ences in colour and metallicity are significant only between the 
extremes of the GC populations. This can be immediately vi- 
sualised using the cumulative distributions of the differences in 
colour and metallicity for the individual stars (i.e., the difference 
between each star and the average of the P population). Fig. [3] 
shows these cumulative distributions, separated for P, I, and E 
stars. The P and I distributions look very similar in both pan- 
els, while the E stars show a distinct behaviour (the probability 
that they are taken from the same parent distribution of the P 
stars are very low: 4x10""* and 0.015 for colour and metallicity, 
respectively): as a whole, they are bluer and metal-richer. 

3. He from the luminosity function of Na-poor and 
Na-rich stars 

With the intent of finding an additional evidence supporting our 
scenario for the correlation between Na-rich stars and their ini- 
tial helium content, we tested also the possibility to detect any 
brightness difference between the RGB bumps corresponding to 
the two distinct stellar populations. In fact, canonical stellar evo- 
lution predicts that the RGB bump brightness i ncreases when the 
initial He content increases (see, for instance 'Salari s & Cassisi 
|2005). As a consequence, one should expect to observe this oc- 
currence when comparing the prim ordial stellar compo nent with 
the intermediate/extreme one (e.g. JSalaris et alil2006l) . 

However, the possibility to detect this observational signa- 
ture in the luminosity function (LF) of a given GC is generally 
strongly hampered by the inadequate statistics (the number of 
RGB stars in the various LF bins is not large enough). We made 
a first explorative attempt (Paper IV) combining the samples 
in NGC 6218 and NGC 6752 to improve the statistics. Here, 
in order to better overcome this problem, we decided to com- 
bine the data for the LFs corresponding to the various GCs in 
our database. We took the values for Ybump from literature (see 
Table nil but could apply the test only to 14 of the 19 GCs. We 
did not reach the level of the bump in some of our target clusters, 
either because they were very far, like NGC 6388, NGC 6441 
or because they were so rich, like NGC 104, NGC 2808, that 
we had enough target in the upper RGB and did not need to 
reach down along the RGB (our original goal was to obtain good 
spectra to measure Na and O abundances for about 100 stars per 
cluster, so we selected preferentially bright targets). Finally, we 
did not find a Yhump value for NGC 7099; since it has the same 
metallicity of NGC 7078, we could have adopted the value for 
the latter, taking into account the distance moduli, but we would 
have introduced an additional uncertainty, given also the slight 
difference in age (Carretta et al. 2009 ^ note th at these two GCs 
have the same age iu iMarin-Franch et al.ll2009h . 
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Fig. 4. Upper panel: Histogram of the difference in magnitude 
between individual V and Vtump (for each cluster) for all stars 
with Na measured (open histogram) and for first-generation (red 
histogram) or second-generation (blue histogram). Lower pan- 
els: Zoom near the RGB bump for all, second generation (I and 
E), and first-generation stars, with Gaussians fitting each bump. 
In each panel the peak of the bump and the associated error are 
indicated. 



We have in total 1368 stars in the 14 GCs with a measured 
Na abundance. To s eparate P stars from I and E ones we used 
[Na/Fe],„,„ -i- 0.3 (see lCarretta et alil2009al) : we have 438 P stars 
and 930 I, E ones (we sum the second-generation populations, 
since the E stars are a minority). We computed for each clus- 
ter the residuals between the individual V values and Vbmtp, to 
eliminate the dependence on metallicity and cluster age. We pro- 
duced histograms both for the whole sample, and for the P and 
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P and I/E components with the same chemical mixture 
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Fig. 5. The V-band RGB bump brightness difference between the 
P sub-population and the I/E one as a function of the He abun- 
dance difference between the two sub-populations. In this case, 
for both sub-populations it has been assumed the same normal, 
ff-enhanced heavy elements distribution (case a in Sect.|2]l. The 
solid line represents the measured empirical RGB bump lumi- 
nosity difference between the P and I/E sub-populations, while 
the long-dashed lines shown the same value at +lcr level. 
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Fig. 6. As in Fig|5] but taking into account the fact that the I/E 
sub-population is characterised by a peculiar chemical pattern 
with a CNONa anticorrelation (case b). 



I/E stars separately. The upper panel of Fig. |4] shows the result- 
ing luminosity functions. Since we have here only the stars for 
which we obtained a spectrum and not the entire RGB, the part 
below the bump is strongly incomplete and we do not see the 
usual increase going towards fainter stars found in photomet- 
ric works on the RGB luminosity function. However, the bump 
is visible both in the whole sample and in the two sub-samples 
which show slightly different peak values. The small offset from 
zero in the histograms is due to the the different photometries 
adopted in the definition of Vbump and in our programme (see 
below). 
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To better appreciate this, we present in the lower panels of 
Fig- Sa zoom in the bump region, separating the three compo- 
nents. Fitting a Gaussian to the bumps we derive a difference in 
magnitude between P and I, E stars of 0.044 mag, in the direction 
expected from theory, with I, E stars presenting a brighter RGB 
bump than P stars. Unfortunately, the significance is low, be- 
cause of the errors on the position of the peaks. A first source of 
error is the statistics: the bump peak (i.e., above the background) 
is produced by only about 150 stars, two thirds of which in the I, 
E population and one third in the P one. To halve the errors in the 
position of the two separate peaks we should increase the num- 
ber of stars by at least 3 times and this is not an immediate possi- 
bility because of the large amount of telescope time and analysis 
implied. Other sources of errors come from the photometry: i) 
Vbump is measured at best with an uncertainty of 0.05 mag, ii) 
the magnitudes we use have a small attached error, and iii) we 
did use a different photometry from the one on which Vhump was 
determined (although the possible offsets should be amply ac- 
commodated within the uncertainty on Vbump)- However, even if 
we tried to re-determine the bump positions in our photomet- 
ric data, the situation would not become significantly better We 
note in f act that the bump ha s an intrinsic width of the order of 
0.2 mag JSalaris et al.L 12002 ) and that we have akeady reached 
that level: the FWHM of the Gaussian used to fit the P population 
(the one with a single Y value) is 0.21 mag (while the FWHM 
of the second-generation stars, with varying Y content, is larger, 
about 0.28 mag). 

In order to estimate the average difference in the initial He 
content between the P stellar population and the I/E one, from 
the measured RGB bump brightness difference, we used the ex- 
tended set of low-mass stellar models for various assumptions 
about the initial He content described in Sect. |2l 

As a first step, it has been assumed that the P and I/E stel- 
lar populations have exactly the same heavy elements distribu- 
tion (case a); then for a fixed age of ~ 12.5 Gyr and iron con- 
tent [Fe/H]=-1.31, we have computed the synthetic LlPI for 
the P population by using a canonical He abundance (Y=0.248), 
whereas for the I/E population various synthetic LFs have been 
computed by using various He-enhanced abundances in the 
range from 0.26 to 0.40. By using these theoretical LFs, we have 
estimated the brightness difference between the RGB bump of 
the P population and that of the I/E one as a function of the dif- 
ference in the initial He contents of the two populations. The 
results of this analysis are shown in Fig.|5] The comparison be- 
tween the model predictions and the estimated empirical differ- 
ence suggests a difference in the mean He content between the P 
and the I/E sub-populations of about 0.01 + 0.01. 

Case b should represent a more significant comparison be- 
tween theory and observations; the adopted heavy elements dis- 
tribution corresponds to a mixture where the sum (C-hN-hO) is 
enhanced by about a factor of 2 with respect to the reference, 
a-enhanced mixture. This value is consistent - although it rep- 
resents an upper limit - with the resul t s of t he spectroscopi- 
cal analysis performed bv^ Carretta et alj (l2005h for the extreme 
values of the chemical anti-correlations observed in GCs. The 
RGB bump brightness difference between the P sub-population 
and the I/E one, when for the latter a CNONa peculiar mixture 
is assumed, is shown in Fig. |6] as a function of the He abun- 
dance difference. This figure deserves some comments: the RGB 
bump brightness of the I/E sub-population is fainter than that 



'" In order to take into account the observational errors that can affect 
the photometric data for each GC in our sample, when computing the 
synthetic LFs we have considered a photometric error of ~ 0.02 mag. 



of the P sub-population for an He content difference lower than 
about 0.045. This occuiTence is due to the fact that when com- 
paring ff-enhanced stellar models with those accounting for a 
CNONa peculiar pattern, we are considering the same iron con- 
tent [Fe/H], and this means that the global metallicity of the stel- 
lar models with the peculiar heavy elements distribution is larger 
than that of the reference models. Therefore, as discussed by 
Pietrinferni et al. (2009) the RGB bump - at fixed age and [Fe/H] 
value - is expected to be fainter in these sub-popul ations char- 
acterised by a peculiar chemical pattern (see also iSalaris et aI.L 
2006). However, for increasing difference in the He content, the 
effect of the He abundance on the RGB bump brightness over- 
comes that associated to the larger metallicity, and the bump of 
the I/E sub-population becomes brighter with respect the P sub- 
population. Theoretical models shown in Fig. |6] suggest an av- 
erage difference between the P and the I/E sub-population of 
~ 0.06 + 0.015. This is in better agreement with the estimates 
obtained in previous sections from independent indicators, even 
if lower. Consider however that NGC 2808 is not in the sample 
and that this cluster brings a strong signature of Y enhancement. 

As stated above, all the stellar models used in the present 
analysis have been computed by assuming a constant [Fe/H] 
value regardless of the adopted initial He content. It is in fact 
expected that both candidate polluters (AGB and FRMS) do not 
alter the initial iron content. However, it could be worthwhile to 
check at what extent - if any - the present results are affected 
by this assumption. For this purpose, we have computed suitable 
stellar models for selected initial He contents, keeping the iron 
content fixed. We found that the magnitude difference between 
the RGB bump of the I/E sub-population and that of the P one 
is increased by a negligible 0.01 mag for a He difference be- 
tween the two sub-populations equal to AY = 0.05 and of about 
0.09 mag for AY = 0.10. From Figs.|5]and|6] it is clear that this 
change does not affect our previous estimates of the He content 
difference. 

Finally, we wish to draw attention on a possible test; clusters 
with larger AY should also display broader RGB bumps. From 
Figs.|5]|6]we estimate that a AY=0.1 should correspond to a vis- 
ible broadening of the bump (since first and second-generation 
stars should display differences of about 0.20-0.35 mag between 
their respective bumps). Precise photometry of all their RGB 
stars should be used. An ideal couple for this comparison are the 
two GCs NGC 2808 and NGC 6121 who have similar metalhcity 
but very different extension of the Na-O anticorrelati on (Carretta 
et al, paper VII) and very different implied AY dGratton et aI.L 
I2010h . 

4. Discussion and summary 

Recently, spreads and splits in GC RGBs other than the ones 
(in V - K, metallicity and temperature) considered here have 
been presented. They are most probably due to the same phe- 
nomenon we are seeing here, since the He-enhanced stars should 
also show other chemical signatures. For instance, the coiTe- 
lated N-enhancement should be responsible for the effects seen 
in the colour-magnitude diagrams involving the Johnson U or 
Stromgren u filters (Marino et al., 2008; Yong et al. , 2008). On 
the other hand, the e ffect seen i n hk p hotometry, attributed to a 
spread i n calcium by Lee et al.l (l2009h . and demonstrated not to 
be so bv lCarretta et al.l(l2010h . still awaits a definite explanation. 
We have used information on RGB stars in GCs to infer the 
plausible He differences implied by the existence of two gener- 
ations of stars. This is a new approach, since the different meth- 
ods used in the past to deduce the He content of GCs gener- 
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Table 4. Adopted values of Vhump- 



GC 



'hump 



Reference 



NGC 0288 


15.45 


0.05 


1 


NGC 1904 


16.00 


0.04 


2 


NGC 3201 


14.55 


0.05 




NGC 4590 


15.15 


0.05 




NGC 5904 


15.00 


0.05 




NGC 6121 


13.40 


0.10 




NGC 6171 


15.85 


0.05 




NGC 6218 


14.60 


0.07 




NGC 6254 


14.65 


0.05 




NGC 6397 


12.60 


0.10 


3 


NGC 6752 


13.65 


0.05 




NGC 6809 


14.15 


0.05 




NGC 6838 


14.80 


0.15 




NGC 7078 


15.25 


0.05 







References. (1) iFerraro et alj l ll999t) : (2) IZoccali et alj ( 119991) : (3) 
lAlves & Saraiedinil(ll99Sf) 



ally involve the HB, the evolutionary phase more sensitive to 
even small He variations (the two exceptions being oj Cen and 
NGC 2808 with their multiple MSs). 

Apart from the information coming from photometric data. 
He can be deduced from spectra of HB stars, where He 
lines are visible for temperatures hotter than about 8500 K. 
I Villanova etsl\ (l2009l) studied seven HB stars in NGC 6752, us- 
ing UVES spectra at high-resolution and very high S/N. For the 
four stars where they were able to measure the He line, they 
obtained a value similar to the cosmological one (Y=0.245). 
However, these four are all O-rich, Na-poor, while Villanova 
et al. could not measure He for the only Na-rich, O-poor star, 
the one that should also be He-enhanced, so the case is not 
yet settled. However, the photospheric abundances of HB stars 
hotter than about 11500 K are altered by atomic diffusion pro- 
cesses and cannot be used directly to measure the original val- 
ues. For instance, Behr (2003) found a depletion in moderately 
hot HB star s in several GCs. M oehler et al., in a series of 
papers (e.g., iMoehler et al.L |2007|) concentrated instead on the 
hottest part of the HB (and blue-hook stars in the two GCs - 
a» Cen and NGC 2808- where they are present) and found He- 
enhancement. This can be attributed either to pollution from a 
previous generation or to He flash induced mixing occuiTing in 
hot He-flashers, i.e. those stars experiencing the He flash not at 
the RGB tip but along the W hite Dwarf cooling sequence (see 
ICastellani & Castellanilll993h. A further di s cussio n of this prob- 
lem in cl) Cen can be found in lCassisi e t al.' ('200j|). 

We have studied in a companion paper ( Gratton et al.L 1201 Ol) 
the effect of (even a small) difference in He content on the mor- 
phology of HBs, and found that He is most probably the third 
parameter governing the HB, after metallicity and age. We re- 
fer to that paper for a long and detailed discussion on the effects 
of He vari ations, on m ethods to measure He content (e.g., the 
R-method: llbenl [l 9681) . on comparison between results obtained 
through different approaches. In the present paper we limit our 
analysis to the differences in He that can be deduced from RGB 
stars. 

May we use the present analysis, combined with the previ- 
ous work on the Na-O and other (anti)correlations to try discrim- 
inating between AGB and FRMS as the polluters of the second- 
generation stars? We need to consider that He is produced es- 
sentially in MS for both classes of polluters, while there is a 
difference for the other "peculiar" elements (O, Na, Mg, Al, Si), 



produced during Hot Bottom Burning or in MS, respectively. 
However, we still miss important details of stellar evolution (e.g., 
how much He is dredged up in AGB stars, or what is the depen- 
dence on metallicity -via rotation- of FRMS), the numbers are 
still small and the inferences are not conclusive. We need to fur- 
ther investigate the dependencies of variations of He and those 
other elements on the stars and cluster properties. 

In summary, with the present work we have seen that: 

i) it is possible to deduce variations of He also from RGB stars 
in GCs when using large samples, treated homogeneously; 

ii) these variations are generally measurable from colour (here 
V - K) and [Fe/H] only between first-generation (P) and ex- 
treme second-generations (E) stars; 

iii) the variations implied for the average of the 19 GCs consid- 
ered are of the order of 0.05-0.10 (AY); however, the heavy- 
elements mixture assumed for the second-generation stars 
has to be taken into account (and we have presented two ex- 
treme cases); 

iv) NGC 2808 is a notable exception: in this cluster, differences 
in colour and metallicity are seen also between P and the 
intermediate (I) component of second-generation stars. This 
is perfectly in line with the fact that NGC 2808 presents a 
clear evidence of three He levels also in MS. For NGC 2808 
the deduced AY values are higher (see Table [3]) than for the 
average of all GCs; 
v) similar results can be obtained also considering the luminos- 
ity of the RGB bump. We could test this method only for 
14 GCs; unluckily, NGC 2808 is not among them. Within 
the limits of the rather poor statistics, we found that the 
implied AY values are in reasonable agreement with those 
found from colour and metallicity; 

vi) the absolute calibration of the AY in GCs is still to be de- 
fined. 
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